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We present preliminary results for the charm quark mass and the D and Ds mesons de- 
cay constants fo and fo, from a lattice QCD calculation with Nf = 2 dynamical fermions. 
We use the twisted mass fermionic action defined at maximal twist so that physical quanti- 
ties are automatically ^{a) improved. Two lattice spacings are considered. The charm quark 
mass has been renormalised in the RI-MOM scheme. After a matching to the MS scheme, we 
obtain from the simulation at a fine lattice {a - 0.09 fm) mp(m,.) = 1.481 ± 0.022 ± 0.092 
GeV, /d = 205 ± 13 ± 17 MeV, /d^ 271 ± 6 ± 6 MeV and from the simulation at the finer 
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larly at the finer one, the error on m^ is dominated by present uncertainty on the renormalisation 
constant Z/>, which should be reduced before performing a reliable continuum limit. 
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1. Introduction 

The physics of charm bound states regained recently the attention of particle physicists with 
the discovery of the new resonances X(3872), X(3943), 7(3940), 7(4260) and Z(3930) [|I|] and 
of a very narrow scalar state 0^(2317) [Q], whose composition is still an open question. The 
experimental evidence for oscillations in the D° — system [^] might be the first signal for physics 
beyond the Standard Model in the charm sector. Even if it is expected that the long-distance 
physics is a dominant effect in that process, the AC = 2 contribution to xd = AMd/Fd might be not 
negligible. It is given by a box diagram, as in the B — B system, and it is proportional to In the 
CKM matrix, Vcs is one of the elements having the largest uncertainty when one does not impose 
the 3x3 unitarity: = 9.82% [Q]: most of it comes from the theory. An appropriate way to 
extract is to measure the leptonic decay width — > Iv, which however requires the estimation 
of the decay constant fo^ . Eventually it is well established that Vd, can be constrained by analysing 
the inclusive semileptonic decay B — > Xclv. An OPE is used in the Heavy Quark Expansion (HQE) 

and the total branching ratio reads ^ BRsl(S) = BR^l^ [1 + /(r)] , r = (^!^^ , where BR^l' the 
dominant term in the mc/mt, expansion. 

From those phenomenological aspects we conclude that it is important to determine as precisely 
as possible the decay constants fo, fn, and the charm quark mass m^. That quark is of course too 
heavy to make predictions by using Chiral Perturbation Theory (;tPT) and too light to compute the 
amplitudes by using only the HQE: corrections of ^(Agco/w^c)" and 0'{m^-/m},)" might be larger 
than the precision of few % that we want to reach. 

Lattice QCD is a good tool to study the charm sector. However the continuum limit can be difficult 
to reach because of large cut off effects (typically, 0.2 < am[! < 0.4). It is therefore crucial to 
improve the action and the currents regularised on the lattice. Several theories proposed in the 
literature so far have the common property that they require the tuning of a certain number of 
parameters to achieve improvement, by applying the Symanzik's program [^. On the other side 
it has been shown that Twisted mass QCD (TmQCD) ||8|, |^ requires the tuning of a single 
parameter (the untwisted mass /mq), so that hadronic quantities like the pseudoscalar meson masses 
and decay constants are automatically (a) improved at maximal twist. Other nice properties of 
such an action are that the physical quark mass is related to the twisted mass parameter of the action, 
its renormalisation is only multiplicative and the pseudoscalar decay constant does not require the 
introduction of any renormalisation constant : for two quark flavors 1 and 2 (for example a light 
flavor i and a heavy flavor h) it is simply given by 

/f5(Mi,M2)= f,^^' . KO|P'-(0)|P)|, P' = %{r)^W2{-r), (1.1) 

where r is the Wilson parameter and we define the composite operator in the physical basis. 

We present a preliminary lattice QCD determination of the charm quark mass nic and the decay 
constants fj) and fj)^. We have performed full dynamical simulations for Nf = 2 light degenerate 
sea quarks; the strange and the charm quarks have been added in the valence sector. The calculation 



is based on the analysis of the gauge configurations ensembles Bi - B5 and C\ -C4 [10| (240 and 
130 configurations of Bi - B5 and Ci - C4, respectively, have been analysed) which have been 
generated with the TlSym gauge action at j3 = 3.9 {a = 0.0855(5)(31) fm) and j8 = 4.05 (a = 
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Figure 1: (a) Chiral extrapolation of the pseudoscalar heavy-light meson masses (j3 = 3.9). (b) Pseu- 
doscalar heavy-light meson mass in function of the heavy quark mass (/3 = 3.9). 



0.0667(5) (24) fm) respectively and the twisted mass fermionic action defined at maximal twist. 
The light quark masses are in the range [ms/6, 2m.s/3], to perform the chiral extrapolation, and we 
used masses around the strange mass and the charm mass to do the appropriate interpolations. 
At each sea quark mass we have computed the two-point correlation functions of pseudoscalar 
mesons. Each measurement has been separated by 20 HMC trajectories, which is enough to avoid 
autocorrelation time effects. The statistical accuracy has been improved by using all to all stochastic 
propagators. Statistical errors on the meson masses and decay constants are evaluated at a given sea 
quark mass by using a jacknife procedure with 10 measurements discarded in each bin. The error 
obtained after a combination of data coming from simulations with different sea quark masses (i.e. 
statistically independent) is computed using a bootstrap method. 



2. Charm quark mass 

To estimate the charm quark mass we use the following strategy: we compute the pseudoscalar 
meson mass mps(/isea)Mf )M/i) (where £ and h are valence light and heavy quarks respectively) at the 
points He = Hsea> which allows us to extrapolate down to the physical light quark mass Hud = i^n^^ 



2 



previously determined in []11[]. Once the dependence on the light quark mass has been taken into 
account, one studies the dependence of mps on the heavy mass /i/j. The bare charm quark mass He 
is determined by using the following condition: mps{Hud-,i^c) = ^d- 

We have performed a quadratic extrapolation of mps{Ht = MseaiM/i) in He down to He = Murf- 
This is illustrated in Figure [T] (left). We have also introduced a logarithmic dependence on Hf- 
fnps{He,Hh) = co{Hh) + [ci(iU/i) +C2(/i/i)ln(a/i£)] . As a third possibility we have done a sim- 
ple Unear extrapolation. The spread between these different fits is included in the systematic error 
at the end of the computation. We then performed a fit of mps as a function of Hh- we used 
either a quadratic polynomial in Hh, or a quadratic polynomial in or also, as a third ansatz, 
mps = do+diHh + We show in Figure [l] (right) the quality of the latter fit function, that appears 
to be the most appropriate to describe the data. Once the bare charm quark mass He is extracted, we 
renormalise it in the RI-MOM scheme: m^-^^^ = Zf-^^^Hc where zf;^-^^^ = \/Zf^-^^^ 
in TmQCD [^]. Finally we perform a matching onto the MS scheme. We indicate in Table |l| the 
value of mp(mc) for the two lattice spacings, knowing that Zf-^°^{I5 = 3.9, 1 /a) = 0.39(1)(2) 
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Figure 2: (a) Chiral extrapolation of mps{lJ.se3, fie — ms,IJ.h ^ nic) (j3 = 3.9). (b) Chiral extrapolation of 
mps{lJ.seiL,Hh,tlh) (J3 =3.9). 



and the preliminary value of Zf^"'^°'^(j8 = 4.05, =0.40(l)(4)i The first error on m,. is 
statistical, the second is the systematic error coming from Zp, the third comes from the uncertainty 
on a and the last one is the systematic error from the chiral extrapolation. In that table we have 
also collected the result of two other determinations of m^, by using the following renormalisation 
conditions: 

1) mp5(/isea = l^ud,l^{ = l^s,l^c) = mo, and 2) mps{n^e^ = jXud , l^c, ^i-c) = m^,- 

Us is the bare strange quark mass which has been determined in [|T4|]. The heavy-heavy pseu- 

doscalar meson correlator has been computed by using the interpolating field i///j(r)7^i///,(r). 

The dependence of those two observables on the sea quark mass is very weak, as shown in Figure ^. 

We used a linear fit in Hi and /^sea to interpolate to /i, and to extrapolate down to H^d respectively: 

'WPs(i"sea,Mf,At/i) = PoilJ.h)+alJ.ePl{lJ.h)+alJ.sea\P2ilJ'h)+aHfP3{Hh)], 

The uncertainty coming from the chiral extrapolation is quite reduced compared to mp5(/Xsea = 
However we introduce a small uncertainty from on the first observable (third error on 
m^^(mc,wJoJ) and from the disconnected diagram which contributes to the second but that we did 
not compute. The last error on m^^{mc,rnD,) and {nic ,m-q^) comes from the uncertainty on a. 
We have not performed a continuum limit extrapolation yet because only 2 lattice spacings have 
been considered so far and the uncertainty on Zp is still rather large, especially at j8 = 4.05. 
Concerning cut off effects, it is remarkable that they appear to be rather weak on the unrenor- 
malised charm mass at j3 = 4.05: indeed, the value of jXc extracted from the 3 observables are 
much closer at this jS than at j3 = 3.9. 

To conclude this section we note that our values of nic are large with respect to most of the recent 
lattice estimations 1 15]: however the currently large uncertainty on Zp at ^ = 4.05 implies that any 



conclusion about the continuum limit result would be untimely. 



'At this /3 a chiral extrapolation has been performed in the valence sector, at a/ijea = 0.003, but not yet in the sea 
sector. However it was found at jS = 3.9 that Zp depends only weakly on the sea quark mass. Thus as a first step of the 
analysis we will include the sea effects in the systematic uncertainty. Moreover, an alternative estimate of Zp at P = 4.05 
may come from scaling as described in jT^], which brings our final estimate of the systematic error to 0.04. 
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Figure 3: Chiral extrapolation of fpsi^i = jUsea, M/i ^ Mc) at )3 = 3.9 (left) and j3 = 4.05 (right). 



3. Heavy-light meson decay constants 

To determine fj) and /^^ we employ the same strategy as in the previous section, using 
eq. (|0|). We show in Figure ^ the chiral extrapolation of fpsii^-xa = i^hp^h) at /X/, ~ /i^ down 
to the physical light quark mass. We found that introducing a quadratic term in the extrapolation 
improves the fit, particularly for the coarse lattice. As before, we also introduced also a logarithmic 
dependence on }X( in the fit. We include the spread between the different chiral extrapolations in 
the systematic error. 

At Hf Us, fps{l^se!L,IJ-e,IJ-h) has a similar Unear dependence on He to the one of nvs(Msea)Mf 
Moreover for both the dependence on /igea is weak and very well described by a linear fit as well. 
We performed a fit of fps\/nips with a quadratic polynomial in ^ (we can not isolate any log- 
arithmic dependence on as{nips) in our range of heavy masses): fps^/mps = /o + + ^ • We 
show in Figure^ the quality of the fit for fps y/mps{ /i„ j ,/!/,). 

We give in Table ^ our values of fo, fo, and foj/o- The first error is a statistical error, the second 
error on fo, and on foJ/o comes from the uncertainty on the bare strange quark mass, the second 
error on fo and the third eiTor on fj)^ come from the uncertainty on the lattice spacing. The last 
error on fo and foJ/o comes from the spread between different chiral fits. We have collected 
the most recent lattice estimations of those quantities in Figure || On the experimental side 
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m^^{mc,mD) 




mf^{mc,mnJ 


3.9 
4.05 


1.481(22)(63)(8)(27) GeV 
1.474(41)(129)(15)(5) GeV 


1.450(12)(61)(10)(15) GeV 
1.498(6)(130)(12)(18) GeV 


1.420(5)(60)(6) GeV 
1.479(2)(129)(8) GeV 



Table 1: Charm quark mass fixed by using different observables. 





3.9 


4.05 


Id 

fn 


205(13)(3)(17) MeV 
271(6)(4)(5) MeV 
1.35(4)(1)(7) 


230(3 1)(6)(6) MeV 
264(5)(4)(7) 
1.13(28)(2)(2) 



Table 2: Decay constants fo and fo, and foj fo from our simulation 
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Figure 4: Scaling law of fpss/mpl{}iudi\^h) a function of l/mps{^ud:^^h) at j3 = 3.9. 
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Alpha, Nf = 0, 2003 
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QCDSF, Nf = 0, 2007 
CP-PACS, Nf = 2, 2000 
MILC, Nf = 2, 2002 
ETMC, Nf = 2, 2007 
FNAL/MILC, Nf = 2 + 1, 2005 
HPQCD/UKQCD, Nf = 2 + 1, 2007 
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Figure 5: Recent lattice computations and experimental measurements of the decay constants fo (left 
panel) and fo^ (right panel). The ETMC values that we indicate are obtained by doing an average of the data 
atj3 =3.9 and J3 =4.05. 



CLEO-c measured fo = 222.6 ± 16.7+^^ MeV Q, /c, = 274 ± 13 ±7 MeV [|T8|], which is a 
combination of analysis of D.v jJ. and Ds T leptonic decays. Note that BABAR measured 
fo, = 283 ± 17 ±7 ± 14 MeV 



4. Summary 

We have presented preliminary results of a lattice computation of the charm quark mass rric 
and the D and mesons decay constants fjj and fo^ by using the Nf = 2 TmQCD action defined 
at maximal twist. Encouraging results are found concerning cut off effects. Indeed, the bare charm 
quark mass extracted from 3 different observables looks consistent at the finer lattice. However, 
before performing the continuum limit on m^, fo and fo^, we still have to reduce as much as 
possible the uncertainty on the renormalisation constant Zp and to increase the statistics at j8 = 4.05. 
A more detailed study of the light quark dependence of the D mesons masses and decay constants. 
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based on heavy-light chiral perturbation theory, is still missing at this stage. Finally a better control 
on the continuum hmit extrapolation will come from the on-going analysis of the data at a coarser 
lattice (a ~ 0.1 fm). 
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